Fourfold degenerate columnar-dimer ground state in square lattice antiferromagnets 
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We construct and study two frustrated quantum spin-1/2 models on square lattice, which are 
like the antiferromagnetic J1-J2 model with some additional four-spin exchange interactions. These 
models admit an exactly solvable case in which the ground state consists of four degenerate columnar- 
dimer singlet (CDS) configurations. Away from the exact case, we employ bond-operator mean- 
field theory to investigate the evolution of the ground state by varying the interaction parameters. 
The mean-field calculation reveals a quantum phase diagram in which the CDS phase undergoes a 
continuous phase transition to either Neel or collinear ordered antiferromagnetic phases. 

PACS numbers: 75.10.Jm, 75.30.Kz, 75.50.Ee, 75.40.Mg 



I. INTRODUCTION 

The low-dimensional quantum antiferromagnets are 
exciting physical systems in which the quantum mechan- 
ics, combined with reduced spatial dimensionality and 
frustration, manifests in a variety of novel properties such 
as the zero magnetic moment and the loss of long range 
order in a spin-liquid ground state, the fractionalized ele- 
mentary excitations, namely spinons, or the spontaneous 
dimerization of spins in the valence bond crystal (VBC) 
ground states 1,2 . Over the last two decades, the ques- 
tion of spin-liquid states and fractionalized excitations 
in two dimensions (2D) has been greatly motivated by 
the high-T c superconductivity in the cuprates. Start- 
ing with Anderson's proposal of the resonating valence 
bond theory of superconductivity 3 , a large body of theo- 
retical investigations is centered around the notion that 
the superconducting state in these materials is derived 
from a spin-liquid state by doping it with charge carri- 
ers. This novel idea continues to inspire the search for 
new quantum spin models with novel ground states on 
two dimensional lattices, and newer methods for study- 
ing these models. 

Even prior to the motivations from the cuprates, there 
has always been a great interest in the studies of quan- 
tum spin models obviously for the reasons of explain- 
ing magnetic phenomena in real materials. Purely theo- 
retical motivations have also led to interesting develop- 
ments in this field. One such notable case is the sponta- 
neous dimerization in the exact, doubly degenerate sin- 
glet ground state of the Majumdar- Ghosh model 4 . It 
is a one-dimensional quantum spin-1/2 model with an- 
tiferromagnetic exchange interactions between the first 
neighbors (Ji) and the second neighbors (J2 = Ji/2). 
This is the simplest known case of a VBC ground state 
whose wavefunction is a pure direct product of the pair- 
wise singlets (also called the valence bonds, or dimcrs) 
between the nearest-neighbor spins. The term VBC 
refers to a spatially ordered configuration of the va- 
lence bonds. This model has inspired a large number of 
studies on frustrated quantum antiferromagnets, partic- 
ularly towards constructing new spin models with exact 
dimer ground states 5 ' 6,7,8,9 . For example, the Shastry- 



Sutherland model 10 (of which there is a material realiza- 
tion in SrCu2(B03)2 11 ' 12 ), the Klein's construction 13 , 
and the AKLT models 14 ' 15 are some of the more notable 
of these subsequent works. 

Through other very important developments in the 
form of Lieb-Schultz-Mattis (LSM) theorem 16 and the 
Haldane's conjecture 17 , it became possible to make more 
general statements on the nature of the ground state of 
quantum antiferromagnets. The LSM theorem states 
that the ground state of a spin-1/2 antiferromagnetic 
chain is unique with gapless excitations, or else it is (at 
least) doubly degenerate. The Haldane's conjecture, on 
the other hand, suggests a unique ground state with a 
spin-gap for an integer spin chain, thus differentiating 
between the integer and the half- integer spins 18 . Sub- 
sequently, there have been efforts to generalize these re- 
sults for different situations 19 . Of particular interest to 
us is a proposal, based on field theoretic analyses 20,21 , 
suggesting the possibility of a fourfold degenerate quan- 
tum disordered ground state in the spin-1/2 quantum 
antiferromagnets on square lattice. In such a disordered 
state, the local magnetic moment is zero due to quantum 
fluctuations, but the lattice translational and rotational 
symmetries can be broken while forming a degenerate 
valence bond crystal state 21 . 

A popular suggestion for such VBC ground states on 
square lattice is the columnar-dimer singlet states (CDS) 
shown in Fig. 2. Since the ground state of the nearest- 
neighbor Heiscnbcrg antifcrromagnet on square lattice 
has Neel order 22 , the CDS phase is expected to emerge 
only when the ordered state is sufficiently frustrated by 
adding suitable competing exchange interactions. The 
nature of the quantum phase transition from an ordered 
antiferromagnetic (AF) phase to a VBC phase is also 
a subject of great current interest. Recently, an inter- 
esting scenario has been proposed in which an ordered 
AF phase undergoes a direct second order transition to a 
VBC phase, through a common quantum critical point 23 . 
In the Landau-Ginzburg approach, the two phases in the 
ground state can be independently characterized by the 
non-zero values of the respective order parameters, and 
the fluctuations thereof. While undergoing a phase tran- 
sition, it therefore seems very unlikely that the two order 



parameters will conspire to continuously vanish exactly 
at the same point. In the new scenario, however, such 
a quantum phase transition is described in terms of the 
spinons, which are 'deconfined' at the quantum critical 
point from which the AF and the VBC phases are claimed 
to emerge as the interaction parameters vary 24 . 

Evidently it is desirable to find a spin- 1/2 quantum 
antifcrromagnet on square lattice which has a fourfold 
degenerate CDS ground state, and where the scenario of 
deconfined quantum criticality can possibly be realized. 
The antiferromagnetic J1-J2 model on square lattice has 
been studied vigorously by various means for exploring 
these possibilities (as well as for independent reasons). 
There was also suggested earlier a certain model with in- 
teractions up to fifth neighbor spins, for which the four 
CDS states were found to be exact eigenstates, but did 
not form the exact ground state 7 . The J1-J2 model has 
Neel order in the ground state when 3\ is sufficiently 
stronger compared to J2, and collinear order (with or- 
dering wavevectors (tt, 0) or (0, n)) in the opposite limit. 
Similar to the Majumdar- Ghosh chain, the J1-J2 model 
on square lattice may be expected to have a sponta- 
neously dimcrized ground state around J2 ~ 0.5Ji. A 
number of studies have indeed shown a spin gapped 
quantum disordered phase in the intermediate range of 
the coupling, 0.4 < J%j J\ < 0.6. This intermediate 
phase is suggested to possess the columnar-dimer or- 
jgj.25,26,27,28,29,30^ j^vrt the alternate possibilities (such as 
a plaquette phase) have also been put forward 31 ' 32 ' 33,34 . 
Here, the nature of quantum phase transition from the 
ordered AF to dimer phase is still under investigation, 
though there are suggestions of its being a weak first or- 
der type 35,36 ' 37 . 

Motivated by these developments, we have constructed 
new models of spin-1/2 quantum antifcrromagncts on 
square lattice. An important feature of the models pre- 
sented here is the case when the CDS configurations form 
an exact fourfold degenerate ground state. This is a case 
of spontaneous dimerization in the ground state, as the 
model Hamiltonians respect the symmetries (point group 
as well translational) of the underlying square lattice. 
Through these models therefore, we are able to realize 
a part of our motivation exactly. The models have two 
parts, one of which is the J1-J2 model, and the remain- 
ing part consists of certain four-spin interactions. For the 
exact ground state, J2 = Ji/2 and the four-spin interac- 
tions are important (but not dominant). Therefore, these 
models can be viewed as two dimensional analogues of 
the Majumdar-Ghosh chain (in the sense of J2/J1, and 
the spontaneous dimerization of nearest-neighbor spins 
in the ground state of desired degeneracy). Away from 
the exact case, we formulate a mean-field theory in terms 
of the bond-operator bosons, to investigate the domain 
of stability of the CDS phase, and the transition to an 
ordered phase. In the following sections, we first discuss 
the models and the exact CDS ground state. Then, we 
describe the mean-field calculations for the general case, 
and present the quantum phase diagram. Finally, we 
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FIG. 1: Pictures-I and II depict the models Hi and Hi, re- 
spectively. The dashed background denotes the J\-J% model. 
The shaded trianlges denote Vs , and the shaded square (in- 
cluding diagonal lines) stands for S' 4 ' (see Eqs. 4 and 6). 
The shaded square in picture-I surrounded by four shaded 
triangles should be understood as a sum of four terms, each 
of which is a product of the square with one of the triangles. 
This sum over the whole lattice gives Hk,i (Eq. 2). In picture- 
II, a pair of shaded triangles with same orientation denotes 
their product, and together they denote a sum of these four 
pairwise products. This gives Hk,2 on square lattice (Eq. 3). 

conclude with a discussion of results, and a summary. 

II. MODELS 

Here, we present two new quantum spin-1/2 mod- 
els on square lattice. Their scheme of construction in- 
volves identifying suitable block-Hamiltonians, and plac- 
ing these blocks (of spins) appropriately on the lattice 
such that the CDS configurations are stabilized as ground 
state. This has been an effective general approach for 
constructing exactly solvable spin models 5 . The Hamil- 
tonians of these models are: Hi = Hj + Hk,i and 
H2 = Hj + Hk,2, in which Hj is the J1-J2 model whose 
Hamiltonian is written as: 

Hj = J\ [S r • S r+a £ + S r • S r + a j;] + (1) 

r 

J2 ^ ^ [Sr • Sr+ax+ay ~\~ S r ' S r _ a £_|_ a y] 
r 

where r is summed over all the sites of square lattice, and 
a is the lattice constant. We parameterize J\ and J2 as: 
Ji = (1 — QJ and J2 = sucn that J\ + J<i = J, where 
J > and < ( < 1. The Hk part of the models is 
new and important as it stabilizes the fourfold degenerate 
CDS ground state. These are written as: 

H K ,i = P f( r ; £) # (4) (r + <i; 6) (2) 

r 8 

h k ,2 = ^EE ? f( r ; £) v \ ( r + £ s ( 3 ) 

where 5 = (S x ,S y ) with 5 X — ±ax and S y = ±ay. The 
coupling K is taken to be positive, and parameterized 
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together with J as: K — k and J = 1 — ft such that 
J + K = 1 and < ft < 1 . These models are pictorially 
described in Fig. 1. 

Operator Va in Eqs. 2 and 3 is a projection operator 
for a triangular block which annihilates Stot = 1/2 states, 
where S to t denotes the total spin of a block. The block 
operator B^ 4 \ of four spin-l/2s, is not a projector, but 
annihilates Stot = sector, and is positive otherwise. In 
terms of the spin operators, these are defined as: 



Va(r; 5) 



■>r+S. 



3 Mr; <£) 



S (4) (r; 5) = - (S r + Sr+i. + S r+5y + S r+Si+Sy 



(4) 
(5) 
(6) 
(7) 



Expanding Eqs. 4 and 6 results in Eqs. 5 and 7 where 
^a(r; £) = (S r • S T+S!C + S r • S r+<5s( + S T+5x ■ S r+Sy ), and 

hi(v, 5) = (S r • Sr+8 X + S r ' S r +6 y + S r • S r+< 5 x+< 5 y + Sr+<5 X ■ 

Sr+8y + Sr+a^ ■ Sr+^+i,, + S r+(5y ■ S r+<5x+(5y ), are the ex- 
change Hamiltonians of the 'completely-connected' three- 
and four-spin blocks respectively 5 . Clearly, V3 is or 1 

for Stot = 1/2 or 3/2 respectively. The eigenvalue of B^ 
is when the spins of a square plaqucttc form a singlet, 
and it is 1/3 or 1 for Stot = 1 or 2 respectively. The Hk,i 
therefore describes the interaction of a square block with 
its neigbhoring triangles, and Hk,i is an interaction be- 
tween triangles. Recently, a model with fourfold degen- 
erate Shastry-Sutherland ground state on square lattice 

has been constructed using a Va-Va interaction 38 , which 

. 2 2 

has inspired the construction of Hk.2 m the present work. 

Rewriting ffi and if 2 in terms of /13 and hi, by putting 
Eqs. 5 and 7 into Hk,i and Hk,2, gives us: H\ t 2 = 
Ho + Vi, 2 + KL, where ifo, V\ and V2 are defined below 
(and, L is the total number of lattice sites; also, ffi, 2 
reads as ffi or ff 2 , and similarly for V\ 2)- 



Ho 



Vo = 



r+ax 



2K 



4if 



EE ft 3(r; 5) M* + §i £) 



r 5 



9 



EE^(r; + 4 fi) 



(8) 



(9) 



(10) 



r S 



Here, Ji = [(1 - $J + §#] and J 2 = [CJ+f-fiT] are 
the effective Ji and J2 exchange interactions. Even when 
the 'external' J is zero (i.e., k = 1), the 'actual' inter- 
actions, Ji and J2, are non-zero, and interestingly their 
ratio is exactly 1/2. Also, in this case, the typical value 
of the four-spin exchange energy (\{Vi t 2)\/L ~ K) would 
be less than (or comparable to) that of the two-spin ex- 
change (\(Ho)\/L ~ 2K) 39 . In this sense, the four-spin 



[in:: 
t r t i 



* * £ * 



4 <j 4 < j 4 



4 <i 4 ^ * 



FIG. 2: Four columnar-dimer singlet configurations which 
form the exact ground state of the models Hk,i and Hk,i- 
A solid line, connecting a filled circle with an empty one, 
represents a dimer (singlet): -^=\ fj — J,f). We denote two 
ends of a singlet bond by a filled and an empty cricle in order 
to emphasise the antisymmetry of its wavefunction under the 
exchange of two ends (spins). 



exchange is important, but not dominant compared to 
the usual two-spin exchange interaction. 



A. Exact ground state 

An important feature of these models is an exact four- 
fold degenerate columnar-dimer ground state that exists 
for k = 1. The models in this case consist purely of 
Hk part, for which J 2 /Ji = 1/2, akin to the Majumdar- 
Ghosh model. For n < 1, the CDS configurations (see 
Fig. 2) do not form the exact eigenstates of ffi and ff2. 
This case will be studied within a mean-field scheme in 
the next section. Presently, we discuss the exact ground 
state. Consider Hk,i on a square lattice in toroidal ge- 
ometry with even number of sites. In this model, every 
square block interacts with its four neighboring triangles 
(see Fig. 1), expressed as a product of B^ with Va. 

Since the lowest eigenvalue of B^ and Vs. operators is 
zero, a 'trial' wavefunction \ifj) for the ground state of the 
model must satisfy the inequality: < E g < (iP\Hk,i 
where E g is the ground state energy. If we can find a 
for which the upper bound of the inequality equals its 
lower bound (i.e., zero), then \ip) is a ground state with 
Eg = 0. The same is also true for Hk.i- 

To construct such an eigenstate, identify the 'local' 
spin wavefunctions which are annihilated by the block- 
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interaction, S( 4 ) ®Vs.. This can be achieved by forming 

a singlet state on a B< 4 > block, or a doublet {Stot = 1/2) 
on a triangle. A singlet state on a square block can be 
formed by making two dimers. Making one dimcr on a 
triangle leaves the third spin 'free', thus forming a dou- 
blet state. A zero energy eigenstate on the full lattice 
can therefore be constructed in the following way: make 
only those dimers which either sit on a square, or on a 
triangle; try forming a configuration of dimers on the 
square lattice such that every block-interaction term in 
Hk,i is simultaneously satisfied (i.e., becomes zero). We 
follow this procedure, considering all allowed dimers at 
the block level, and find the four columnar- dimer config- 
urations, shown in Fig. 2. These four CDS states clearly 
form the exact ground state of Hk,i, because every block- 
interaction term in Hk,\ annihilates them. Following the 
same procedure, we find that the four CDS states also 
form the exact ground state of Hk,i- In this case, the 
interaction between triangles becomes zero. 

The CDS wavefunctions can explicitly be written as: 

LJ2 Ly 

\CDS Xtl ) = <g> {J J] [(2Z-l,m), (2l,m)] (11) 

1 = 1 m=l 

\CDS y ,x) = d \CDS x ,x) (12) 

\CDS X , 2 ) = Cl\CDS Xil ) (13) 

\CDS y , 2 ) = el\CDS xA ) (14) 

where [(711,7711), (712,7712)] denotes a singlet between two 
spins with coordinates (tt,i,toi) and (712,7712) on square 
lattice (L x = L y , and L = L x L y ) 40 , and C4 denotes 
the rotation by ir/2. The subscripts x or y in the state- 
labels denote whether the dimers arc formed along x or in- 
directions, and 1 or 2 denote two states of x or y type. For 
reference, we can take the bottom-left picture in Fig. 2 as 
\CDS Xt i). Clearly, the CDS states break the rotational 
symmetry, and also the lattice translation, as \CDS Xt \) 
transforms to \CDS X .2) under the translation along in- 
direction and \CDSy_i) transforms to \CDS y> 2) under 
the translation along y-direction. Also, the CDS ground 
state has a gap to the spin excitations. The energy gap to 
a localized triplet excitation, created by breaking a sin- 
glet bond, is estimated to be 8/3 (for n = 1). For re < 1, 
the triplet excitation is expected to hop, and therefore, 
the spin-gap will be reduced. 

In principle, one could envisage the possibility of a 
'fifth' state in the exact ground state. To show rigorously 
that the four CDS configurations are the only states in 
the exact ground state is a hard problem. However, as 
stated above, the explicit construction gives only four 
CDS configurations in the dimcr ground state. We sup- 
port it by the following counting argument 38 . On a lat- 
tice of L sites, there are exactly L/2 dimers in any va- 
lence bond state. Since there are exactly 4L number of 
block-interaction terms (both in Hk,i and Hk,2): there- 
fore on the average, a dimer is to be shared between eight 
block-interaction terms. There are only two types of al- 
lowed dimers: 1) between first-neighbors, and 2) between 



second-neighbors. A second-neighbor dimer is common 
only to four block-interaction terms, whereas a first- 
neighbor dimer is shared exactly between eight block- 
interaction terms. Therefore, a dimer eigenstate of the 
models must be made of nearest-neighbor dimers only. 
Also, it is necessary that only one of the two interacting 
blocks becomes zero, otherwise at least one other block- 
interaction term will remain unsatisfied 41 . Furthermore, 
there are four ways in which a given spin can form a dimcr 
with its nearest- neighbors (two each along x- and y-axes). 
Now, if we construct a dimer ground state observing all 
these constraints, we only generate four CDS configura- 
tions (one each for the four choices of the initial dimer). 
This simple argument also demonstrates the dependence 
of the degeneracy of a dimer ground state on the (rele- 
vant) coordination of a spin 42 . For example, the fourfold 
degeneracy of the Shastry-Sutherland ground state in a 
recent model 38 is consistent with the second-neighbor co- 
ordination on square lattice being four. 



III. TRIPLON MEAN-FIELD THEORY 

Now, we study the general case of H\ and H2 in which 
k and £ take arbitrary values in the interval [0,1]. As 
mentioned earlier, the spin-gapped CDS state is no more 
the exact eigenstate of Hi t 2 for re < 1. However, we 
expect that the spontaneous dimerization in the ground 
state will continue to persist until the spin-gap is closed. 
We therefore want to investigate, within some approx- 
imate but physically acceptable scheme, the domain of 
existence of the CDS phase in the re-£ plane. For this 
purpose, we choose to work with bond-operator repre- 
sentation of spin-1/2 operators 29,43 , which provides an 
appropriate tool for studying a spin system with dimer- 
ization. We formulate a simple mean- field theory of H\_ 2 
in terms of the bond-operators, and calculate the quan- 
tum phase diagram within this approach. An important 
feature of this formulation is that we don't have to guess 
which antiferromagnetically ordered phase will the CDS 
phase go into. It is decided simply by the wavevector 
at which the spin-gap closes, and the triplons condense. 
A triplon is a triplet quasi-particle residing on a bond. 
We will formulate this mean-field theory in the present 
section, and the results of the self-consistent calculation 
will be discussed in the next section. 



A. Bond-Operator Representation 

The Hilbert space of a pair of S = 1/2 spins, say Si 
and S2, is spanned by a singlet state |s) and three triplet 
states \t x ), \t v ) and \t z ). In the bond-operator represen- 
tation, these spin states are created out of a vacuum |0) 
by means of a 'singlet' and three 'triplet' creation oper- 
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ators, as defined below. 



B. Mean-Field Theory 
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:=4|0) 


(18) 



The bond-operators and v a (a = x, y, z) are canonical 
bosons with a constraint: s^s + t^ a t a = 1, where repeated 
Greek index is summed over. In terms of these bosons, 
the spin operators, Si and S2, can be expressed as: 

Sia = ^ (sH a + t\s — ie a p y iptjj (19) 
S 2a = i (-sh a - t^s - iecpyt^t^ (20) 

where e a p 1 is the totally antisymmetric tensor. Eqs. 19 
and 20 satisfy the spin algebra (by making use of the 
boson commutation relations and the constraint), and 
hence define the bond-operator representation 29 of two 
5 = 1/2 spins. And, Si • S a = -f s f s + J4*a- 

In principle, one could just rewrite any spin-1/2 model 
in terms of the bond-operators. But this alone does not 
help, because now we have a constrained interacting bo- 
son problem. However, within a physically meaningful 
approximation, this representation can still be used for 
studying interesting problems where we have a scope for 
dimerized singlet ground states. The physical setting, in 
qualitative terms, of such a problem is simple: we have 
a background of dimers with some mean singlet ampli- 
tude per bond, in which an elementary triplet excitation 
can be created by breaking a singlet bond; this triplet, 
in general, can disperse through the background of sin- 
glets, assisted by the exchange interactions present in 
the system. We call this dispersing &on<i-triplct a triplon 
(in contrast to a magnon, which is also a spin-1 excita- 
tion, but created by a single spin- flip). In terms of the 
bond-operators, this implies a simple mean-field theory, 
in which we replace s and by a c-number s, which 
quantifies mean singlet amplitude on every dimcr. The 
resultant Hamiltonian is a model of interacting triplons, 
which can further be simplified by ignoring the interac- 
tion (similar to the spin-wave analysis of a magnetically 
ordered ground state). This two-step simplification re- 
sults in a problem which can be exactly diagonalized. 
On each dimer locally, the bosons must also satisfy the 
constraint. This is however done, on average, by intro- 
ducing a global chemical potential. Besides the dimerized 
ground states, we can also describe magnetically ordered 
phases in the same formulation, by the condensation of 
triplons. Therefore, in this mean-field approximation, we 
can also study quantum phase transitions in the ground 
state, by varying suitable parameters in the problem. 



We now derive the mean-field triplon Hamiltonian for 
Hi 2- Since there are four CDS configurations in the 
exact ground state, we can choose any one of these as 
a reference state for bond-operators. We take \CDS x ,i) 
as the reference state, and represent the spin models in 
terms of the bond-operators. Though we would like to 
work with all four CDS states at the same time, but we 
can not. This is a limitation of the present approach. It 
will, nevertheless, give us a good idea of what happens 
to the CDS phase by varying n and £. 

We substitute Eqs. 19 and 20 for the spin operators in 
Ho, V\ and V2, and incorporate the constraint on bond- 
operators by including an extra term: (Ur(1 — s r s r — 
t\i a t Ra ), in the Hamiltonian (here, R denotes the posi- 
tion of a dimcr). The Lagrange multipliers are taken to 
be uniform, that is //r = /1. By making the simplifying 
approximations discussed above, and Fourier transform- 
ing the triplon operators, we get the following trilpon 
Hamiltonian for both Hi and H 2 - 



H t = £o + ^{(A-.s 2 £ k ) 



l ka -ka ' t -kQ t kQ 



where 



\ = \{i-QJ+\k-h 

is the effective chemical potential of triplons 
(1 



} (21) 



(22) 



En = - 



- C) T &K 5A x 2 



(i 



e L (23) 



is a constant term, and is written as: 
^k 



J{(1 - C) (cos2fc K - 2cosfc y ) /2 + 

,} + —(!- s 2 ) (cos2fc a 



C (1 + cos 2k x ) cos k y ] 
— cos k y + cos 2k x cos k y ) 



(24) 



Remember that K = k, J = 1 — K, and L is the total 
number of lattice sites. Also, the lattice constant a has 
been put equal to 1. The Fourier transformation of the 



triplon operators is defined as: tR Q = y jr J2 k e lk R tka, 

where k belongs to the CDS Brilluoin zone (CDS-BZ) 
which is folded half-way along x-axis with respect to the 
Brilluoin zone of square lattice (Sqr-BZ) (see Fig. 3). 

We diagonalize H t by making Bogoliubov transforma- 
tion: t ka = 7 kQ cosh6> k +7l kQ sinh# k , where 7 kct are new 
bosons. Rewriting H t in terms of 7 kQ , and demanding 
the off-diagonal terms to be absent gives the equation: 
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Here, ?7k = £k 



^s 2 (cos2k x 



cos ky + cos 2k x cos k y ). 
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FIG. 3: The shaded rectangular region denotes the Brilluoin 
zone of the CDS state (CDS-BZ), while the dashed-outlined 
big square is the Brillouin zone of the underlying square lattice 
(Sqr-BZ). The Neel ordering wavevector Qjv = (7T,7r) in Sqr- 
BZ maps to (0, n) in CDS-BZ, and similarly, the collinear 
ordering wavevector Qc = (tt,0) maps to (0,0). 



tanh 2^k 



2 a/(A 



2 ^k), for 6*k- And, we get the 



following diagonal form of the mean-field Hamiltonian 

3 



(25) 



Here, = \J A(A — 2s 2 £fc) > is the triplon quasi- 
particle dispersion. There are three degenerate quasi- 
particle branches (each corresponding to a different a). 
The vacuum of the 7 kce bosons is the ground state of H t , 
with ground state energy density, e g : 



These self-consistent equations are a saddle point condi- 



tion on e g [X, s ], that is: de g /d\ — and de g /ds 2 



0. 



respectively. As long as the triplon gap is non-zero, 
the denominators in Eqs. 27 and 28 would not have 
any problem, and the two equations can be solved self- 
consistcntly. 

The s 2 thus calculated directly measures the expecta- 
tion value of the singlet projection operator on a dimcr 
in the mean-field CDS ground state. The singlet pro- 
jector on a dimer is defined as: j — Si ■ S2, which in 
terms of the bond-operators is exactly equal to s^s, and 
therefore, has s 2 as its mean-field value. In fact, it is 
encouraging to state that our mean-field calculation is 
exact for k = 1, in the sense that it gives s 2 = 1, which is 
same as the expectation value of the singlet projector in 
the exact CDS ground state. The mean-field theory also 
gives zero magnetic moment in the gapped case, as ex- 
pected for a quantum disordered state. It can be shown 
by noting that (S la ) = s(4+0/ 2 

(*8* 7 ) in the 

mean- field theory, and similarly (S^a)- Since there is no 
mixing of modes with different values of a (see Eqs. 21 
and 25), and the ground state is a 'normal' boson vac- 
uum (i.e., (t a ) = (tl) = 0), the local magnetic moments, 
(Si) = (S2) = 0. Hence, the mean-field ground state 
is quantum disordered, when triplons are gapped. The 
actual dependence of the gap and s 2 on k and £ is com- 
puted in the next section. Below, we discuss the case 
when triplon quasi-particles become gapless. 
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2. Ordered Antiferromagnetic Phase 



where e = E /L. Minimizing e g with respect to A and 
s 2 will give us the mean-field solution. Since we are in- 
terested in understanding the competing phases in the 
ground state, we restrict our investigation to T = 0. 
Though, in principle, we could also formulate a finite 
temperature mean-field theory. 



1. Gapped CDS Phase 

In the CDS phase, it costs finite energy to create a 
triplon. This manifests as a gap in the quasi-particle 
energy, that is E^ > 0. The mean-field parameters s 2 
and A in this ground state are determined by itcrativcly 
solving the following self-consistent equations. 



5 3 ^ A - s 2 £ k 



Eu 



k 



(27) 
(28) 



Suppose, for some values of K and £, the dispersion 
i?k touches zero at a certain wavevector Q. Then, the 
corresponding terms in Eqs. 27 and 28 become singu- 
lar, and there appears a third equation: Eq = 0, relat- 
ing two unknowns s 2 and A in the problem. Following 
Einstein's approach, we interpret the singularity in the 
self-consistent equations as Bosc-Einstcin condensation 
(BEC) of triplons in mode Q. This gives us a 'third' 
unknown in the form of triplon condensate density n c , 
which also resolves the problem of closure of the number 
of unknowns to the number of equations 44 . The physical 
consequence of triplon condensation, as shown later, is 
the emergence of non-zero local magnetic moment, giv- 
ing rise to antiferromagnetic order in the system. 

The closing of the quasi-particle gap fixes A through 
Eq = 0. The singular terms with Eq as denominator 
can be suitably expressed in terms of the condensate den- 
sity, n c . The constraint on the bond-operators implies: 
1 — s 2 . Therefore, the density of con- 
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densed triplons in the ground state can be written as: 



= 7<*Jja*Qa> = 1 -* a - 7 



It. 



a - s 2 & 



E, 



Q 



(29) 



Above, in the last step, Eq. 27 has been invoked. This 
is the rule for associating singular terms to n c . If there 
is more than one Q for the same condensate, it is under- 
stood that we sum over Q in Eq. 29. We now get the 
following self-consistent equations for A, n c and s 2 : 



A 



2s 2 Cq 
1 

5 
2 



k#Q 



1 V — 



s 2 & 



(30) 
(31) 

(32) 



Eq. 30 is a non-trivial solution of Eq = for A (there 
is a trivial but physically inconsitent solution in A = 0). 
Eqs. 31 and 32 are the restatements of Eq. 28 and 29, 
respectively. The scheme of iterative calculation in this 
case goes as follows. For some initial value of s 2 , calculate 
A. Feed this A and s 2 into Eq. 31 to calculate n c . Then, 
calculate new s 2 from Eq. 32, and close the loop. 

In the BEC ground state, the operators iq a , iq Q 
can be treated as a c- number t such that: n c = 
2(^Qc^Qa) I L = where t is taken to be real, and in- 

dependent of a (because three branches are degenerate). 
For the spins of a dimer, we can show that: (Sr.i q ) = 

-(S*r, 2q ) = at cos(Q • R) = cos(Q • R). 

Since the different components of a spin have same ex- 
pectation value, we define the local magnetic moments 
M R1 and M R2 as: (fi^ia) = Ma,i/v3, and similarly 
for M Ri2 . Therefore, Mr.i = -Mr )2 = s^nZ cos(Q-R). 
Having a non-zero n c thus enables the emergence of an 
AF order with an ordering wavevector Q, and an order 
parameter (staggered magnetic moment), M s = s^/n^. 



IV. RESULTS AND DISCUSSION 

We first compute s 2 and A in the gapped phase. For 
k = 1, we get s 2 — 1, which is same as the exact result. 
Also, the flat triplon dispersion with an energy gap of 
2.66 is consistent with our estimate of 8/3 for the gap for 
a localized triplon in the exact case. We further compute 
the mean-field variables for different values of k and £. 
For sufficiently small k, our system undergoes a quantum 
phase transition to ordered AF ground states by closing 
the triplon gap. The locus of those points in the k-£ 
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FIG. 4: The quantum phase diagram of the models Hi and 
Hi in the triplon mean-field theory. The exact CDS ground 
state corresponds to k = 1 line. Also, shown are the regions 
of Neel and collinear ordered AF ground states for £ < 0.16 
and £ > 0.4, respectively. 



plane, where the triplon gap just vanishes, defines the 
phase boundary of the gapped CDS phase. Fig. 4 shows 
the quantum phase diagram of Hi and i?2, as given by 
the triplon mean-field calculation. To correctly identify 
the ordered phases, we must know the minimum of E^ 
in the gapped phase. We find two different values of Q, 
for which the quasi-particle dispersion is minimum (see 
Fig. 5). These are: Q = Q N = (0, tt) and Q c = (0,0) 
(see Fig. 3). The minimum value of the triplon disper- 
sion defines the spin-gap: A = ^/A (A — 2,s 2 £q), in the 
CDS phase. The shift of the minimum of E^ from Qtv 
to Qc occurs precisely at £ = 1/3 (which correspond to 
Jij J\ — J2/J1 = 1/2). The k-( plane has two differ- 
ent regions corresponding to Eq n or Eq c being zero, in 
which the ground state has Neel or collinear order, re- 
spectively. The AF ground states as they emerge on top 
of the CDS background are pictorially shown in Fig 6. 
This figure also shows how Qn and Qc correspond re- 
spectively to Neel and collinear phases. From the gapless 
side, where n c (or M s ) is the relevant order parameter, 
the phase boundaries in Fig. 4 are to be understood as 
those points at which n c continuously vanishes. 

In Fig. 7, we present A and n c as a function of n for 
different values of (. For £ = 0, the gap vanishes at K « 
0.25, and n c continuously builds up for lower values of 
k. This describes a continuous quantum phase transition 
from the CDS to Neel phase. For < ( < 0.16, we always 
find a k, below which the ground state is Neel ordered. 
For 0.16 < C ^ 0.4, the gapped CDS phase persists for 
the entire range of n. For example, the non-vanishing A 
for £ = 0.32 in Fig. 7 shows the absence of a quantum 
phase transition. For £ > 0.4, the CDS phase gives way 
to collinear AF phase for sufficiently small values of n. 
For £ = 1, the transition to the collinear ordered phase 
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(71/2,0) (7t/2,7t) 



(0,71) (71/2,71) (0,0) 



FIG. 5: The dispersion i?k of the triplon quasi-particles. No- 
tice that the minimum of each curve is at Qjv = (0, 7r) for 
C < 1/3 and Q c = (0,0) for £ > 1/3. 

/ / / / / Z / Z / Z 

/ / / / / / / / / / 

Jf / / / / / / / / / 

Jf / / / / / / Jf~J{ / 

/ / / / / / / / / / 

FIG. 6: Schematic drawings of Neel (left) and collinear (right) 
AF ground states in the triplon mean-field theory. Note that 
all dimers have identical magnetic moments in the collinear 
state, hence Qc = (0, 0). In the Neel state, the magnetic mo- 
ments on neighboring dimers along y-direction flip, therefore 
Qjv = (0,tt). 



happens at K ks 0.57. The transition to collinear phase 
in general looks more pronounced as compared to the 
CDS to Neel transition. In Fig. 8, we present the same 
information in terms of the magnetic order parameter, 
M s = Sy/n^. For completeness, in Figs. 9, 10 and 11, 
we also present the data for s 2 , 
range of K for a few values of Q 

We notice that the dimer order of the CDS phase does 
not wash out completely in the AF ground states. That 
is, the chosen dimers always have more singlet ampli- 
tude as compared to other equivalent bonds (pictorially 
represented in Fig. 6). This is a limitation of the bond- 
operator mean-field theory, because only one of the four 
degenerate states can be chosen to carry out the triplon 
analysis. Look at the reduction in s 2 in the ordered AF 
phases with respect to its value in the CDS phase in 
Fig. 9. Clearly, this reduction in the Neel phase is much 
less compared to the collinear phase. In the collinear 
phase at k = 0, the value of s 2 is in fact very close to 
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FIG. 7: The triplon condensate density n c , and the normal- 
ized spin-gap A/A(« = 1) versus s. A/0 gives CDS phase, 
and n c 7^ gives magnetic ordering. 




FIG. 8: The staggered magnetic moment M a and normalized 
spin-gap vs. k, for different values of £. 



a crude estimate of 0.5 (for a classical pair of opposite 
spin-l/2s). This can be qualitatively explained by ob- 
serving that the collinear phase is also columnar like the 
CDS state, whereas the Neel state has a different mag- 
netic lattice. Therefore, undergoing a transition from 
CDS to Neel phase needs some reorganization, which is 
unfortunately restricted in the bond-operator mean-field 
theory. This difference also reflects in the bounds of the 
intermediate phase for k = 0, where Jij 3\ ~ 0.19 (on 
the Neel side) is considerably smaller compared to the 
known numerically esitmates 28 . An alternate approach 
for treating the degenerate dimer states on equal foot- 
ing is very much desirable. That will not only improve 
the results quantitatively, but also help to investigate the 
possibility of spinons in these models. 

The existence of the CDS ground state in H\ and H 2 , 
starting from the exact case and extending all the way 
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FIG. 9: The singlet amplitude on a dimer, s 2 vs. k. 



into the intermediate range 0.16 < C ^ 0.4 (equivalently, 
0.19 < J2/J1 S 0.67) of the pure J1-J2 model, supports 
the view that there is a spontaneously dimerized colum- 
nar phase between the Neel and collincar ordered AF 
ground states in the J1-J2 model. We believe that the 
ground state of pure J2/J1 — 1/2 model is smoothly 
connected to the exact CDS ground state of k = 1 model 
(like a fixed point Hamiltonian). In other words, an ex- 
act numerical study of Hi and H2 is expected to show 
an adiabatic evolution of the exact CDS ground state to 
the actual ground state of pure J2/J1 = 1/2 model. To 
conclude, the mean-field theory presents a qualitatively 
robust quantum phase diagram, which may however be 
improved quantitatively by more exact numerical meth- 
ods. Also, it reveals a continuous quantum phase transi- 
tion between the CDS and AF phases, which should fur- 
ther be ascertained by more exact results. In future, we 
plan to carry out exact numerical studies on these mod- 
els to address many such questions. Finally, the models 
Hi and H2 have the necessary ingredients for discussing 
the possibility of deconfined quantum critical points, and 
therefore should further be investigated in that direction. 



V. SUMMARY 

We have constructed new models of spin-1/2 quantum 
antifcrromagnets on square lattice, with a fourfold de- 
generate columnar-dimer ground state. These models 
consist of a block-interaction part plus the J1-J2 model. 
The block-interactions cither involve two triangles, or a 
triangle and a square. The interactions in the models 
are parametrized as k and C, such that n = 1 gives pure 
block-interactions while k = is pure J1-J2 model, and 
£ is a measure of J2 with respect to J\. For k = 1, the 
ground state is exact, and consists of four CDS configu- 
rations. Away from the exact case, we have formulated a 
triplon mean-field theory, in terms of the bond-operator 
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FIG. 10: Effective chemical potential A vs. k. Change in the 
sign of slope accross the transition suggests a second order 
transition. There is no such change for C, = 0.32. 
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FIG. 11: Ground state energy per site vs. n. 



bosons, to investigate the domain of stability of the CDS 
ground state. The CDS phase is characterized by a gap in 
the triplon dispersion. Since k stabilizes the CDS phase 
while C (close to or 1) favours antiferromagnetic or- 
der, the ground state shows quantum phase transition 
between the ordered AF and the CDS phases in the k-( 
plane. The onset of an ordered phase is characterized by 
the closing of the triplon gap. In the gapless phase, the 
triplons Bose condense, and give rise to an ordered anti- 
ferromagnetic ground state, with an ordering wavevector 
Q given by those points in the Brillouin zone at which 
the triplon dispersion touches zero. 

The mean- field theory is exact for k = 1, as it gives 
a perfect CDS phase. For sufficiently small k, two 
types of ordered AF phases emerge: the Neel and the 
collincar. The Neel phase exists in a bounded region 
close to (k, C) = (0,0), and there is a region of collincar 
order around (k, () = (0, 1). The rest of the phase dia- 
gram is occupied by the stable CDS phase (see Fig. 4). 
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The triplon gap in the CDS phase continuously goes to 
zero while approaching the phase boundary. The same 
is true for staggered magnetization in the ordered AF 
phase. The triplon mean-field theory thus reveals a pic- 
ture of various phases and the transitions thereof, in the 
ground state of these models, which should further be 
tested against and improved by more elaborate numeri- 
cal calculations. 
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